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Abstract 

Different numbers of self-gravitating particles (in different types of periodic motion) are most 
likely to generate very different shapes of gravitational waves, some of which, however, can be 
accidentally almost the same. One such example is a binary and a three-body system for Lagrange's 
solution. To track the evolution of these similar waveforms, we define a chirp mass to the triple 
system. Thereby, we show that the quadrupole waveforms cannot distinguish the sources. It is 
suggested that waveforms with higher £-th multipoles will be important for classification of them 
(with a conjecture of £ < for particles). 

PACS numbers: 04.30.Db, 95.10.Ce, 95.30.Sf, 04.25.Nx 
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Introduction. — Can one see an apple fall at dark night? This is an inverse problem in 
gravitational waves astronomy. It can be specifically stated as "how can we know the 
source information such as the number of particles, their geometrical shape and motion 
from observations of gravitational waves?" This problem is analogous to the well-known one 
for the sound, which was raised by Kac in his celebrated paper |l| entitled "Can one hear 
the shape of a drum?^ Seeking an answer is beyond the scope of this paper. As a specific 
issue which is related with the inverse problem, we shall examine gravitational radiation by 
a certain class of orbital motion of self-gravitating objects. 

In the near future, direct detections of gravitational ripples (and consequently gravita- 
tional waves astronorny) will come true owing to a lot of efforts by the on-going or designed 
detectors mm , la, 13]. One of the most promising astrophysical sources is inspiraling 



re. 
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ativity has succeeded in simulat- 



12( |. Analytic methods also have 



and finally merging binary compact stars. Numerica 
ing merging neutron stars and black holes 
nicely prepared accurate waveform ternplates for inspiraling compact binaries, notably by 
the post-Newtonian approach (See 13l, llJ] for reviews) and also by the black hole perturba- 
tions especially at the linear order in mass ratio (See also [15] for reviews). Bridges between 
;he inspiraling stage and the final merging phase are currently under construction (e.g.. 



16|,ll7|). 



There is a growing interest in potential astrophysical sources of gravitational waves in- 
volving 3-body interactions (e.g., 18, lol and references therein). It is well-known that even 
the classical three-body (or N-body) problem in Newtonian gravity admits an increasing 



number of solutions 
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21| . Some of the orbits are regular, while the others are chaotic. 



For simplicity, we focus on several periodic orbits of three body system; Lagrange's triangle, 
Henon's criss-cross and Moore's figure-eight, which are explained later (See also Fig. [1]). 
Here, it should be noted that Nakamura and Oohara 22| studied numerically the luminosity 
of gravitational radiation by N test particles orbiting around a Schwarzshild black hole, as 
an extension of Detweiler's analysis of the = 1 case by using Teukolsky equation 24], 
in order to show the phase cancellation effect, which had been pointed out by Nakamura 



and Sasaki 



25| . Their N particles are test masses but not self-gravitating. Another inverse 



problem of reconstructing the gravitational wave signal from the noisy data acquired by a 



network of detectors has been discussed (e.g. 
different from those of the existing works. 
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27(1). Our aim and setting are completely 
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The purpose of this paper is (1) to point out a case where very similar shapes of waves 
are generated accidentally by different numbers of particles and (2) to show that the usage 
of higher multipole contributions will be necessary for distinguishing such sources. In order 
to track the evolution of the waveforms, we shall define the chirp mass so as to extend to a 
three-body system. Thereby, we shall show that the octupole order is required to disentangle 
such very similar waveforms that coincide with each other at the quadrupole level. This will 
suggest that theoretical waveforms including sufficiently higher i-th order multipole will be 
important for classification of sources generating such similar waveforms (with a conjecture 
about i and A^). 

Throughout this paper, we take the units of G = c = 1. 
Some periodic orbits for three-body systems. — For simplicity, we assume that the motion of 
massive bodies follows the Newtonian equation of motion. It is impossible to describe all the 
solutions to the three-body problem even for the 1/r potential, as mentioned above. The 
simplest periodic solutions for this problem were discovered by Euler (1765) and by Lagrange 
(1772). The Euler's solution is a collinear solution, in which the masses are coUinear at every 
instant with the same ratios of their distances. The Lagrange's one is an equilateral triangle 
solution in which each mass moves in an ellipse in such a way that the triangle formed by 



the three bodies revolves. Let us take as anot 
orbit found by Henon in 1976 28j (See also 



ler interesting solution the so-called criss-cross 



2^ for the initial condition for each mass and 



recent extensions of the solution) 



Since the figure-eight solution was found first by Moore by topological classification [30[ , 
choreographic solutions have recently attracted increasing interests in astronomy, mathe- 
matics and physics, where a solution is called choreographic if every massive particles move 
periodically in a single closed orbit. The figure-eight solution is that three bodies move pe- 
riodically in a single figure-eight JsO]. The existence of such a figure-eight orbit was proven 
by Chenciner and Montgomery 3l|], where the numerical initial condition for each mass is 
also given. This odd solution is remarkably stable in Newtonian gravity {32, 33[^ Heggie 
discussed a formation mechanism as an outcome from scattering of two binaries ^] . Even- 



tually its unicity up to scaling and rotation has been recently proven 



351]. The trick figure 



eight remains true even if we consider the general relativistic effects at the post-Newtonian 



order 



36l | and also at the second post-Newtonian one [19|]. This is a marked contrast to a 



binary case, which produces a complicated flowerlike pattern by the periastron advance in 
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FIG. 1: Orbital shapes, (a) Top left: Circular orbit for two-body system as a reference, (b) Top 
right: Triangle solution by Lagrange, (c) Bottom right: Criss-cross orbit by Henon. (d) Bottom 
left: Figure-eight trajectory by Moore. 

Einstein gravity. It is interesting to investigate relativistic effects on various kinds of orbital 
motions, which are discussed mostly in Newtonian gravity. It is a topic of future study. The 
radiation by the figure eight has been also investigated [18 1 

Gravitational waves. — In the previous part, we have mentioned several periodic solutions. 
Figure [2] shows the gravitational radiation by massive particles in these periodic motions, 
where the quadrupole formula is used. 
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FIG. 2: Gravitational waveforms in arbitrary units (T = orbital period). Dotted blue and solid 
red curves denote + and x modes, respectively, (a) Top left: Gravitational waveforms by binary 
system with a mass ratio of 2:3 in circular motion, (b) Top right: Lagrange's triangle solution 
for a mass ratio of 1:2:3. (c) Bottom right: Henon's criss-cross, (d) Bottom left: Moore's figure- 
eight. Criss cross and figure eight have larger curvatures in the orbital shapes than Keplerian and 
Lagrangian orbits, which lead to larger acceleration of the particles and thus relatively stronger 
radiation. 

Interestingly, the waveforms from a binary in circular motion and a three-body system 
constituting the Lagrange solution are the same in shape. It is worthwhile to mention that, 
if the third mass is extremely small, its contribution to the quadrupole waves becomes linear 
but not cubic in mass because its orbital radius is of the order of a triangle's side length, 
namely bounded from above. If one adjust properly distance r from an observer to the source 
with the same orbital period, the waveforms (including the amplitudes) could perfectly agree 
with each other. 
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FIG. 3: Definition of Op in the Lagrange's triangle solution. The angle 9p is measured from X-axis 
to the direction of each mass at the initial time. 

Chirp mass for three-body systems. — The waveforms shown above are valid only in short 
term. The gravitational waves will gradually carry away the system's energy and angular 
momentum, and will eventually shrink the orbital size. Consequently, the amplitude and 
frequency of the waves will become larger and higher, respectively, with time. For a binary 
case, the frequency sweep is characterized by its chirp mass. 

Here, we investigate the evolution of the waveforms for a three-body system for the 
Lagrange's solution (on x-y plane). The initial positions of each mass denoted by rrip {p = 
1,2,3) are expressed as a;i = (0,0), X2 = a(v/3/2,l/2), and a^s = a(0, 1), where the side 
of a regular triangle is denoted as a. We take the coordinates such that the center of mass 
(COM) is at rest as {xcomi Ucom) = a(v^z/2/2, (z/2 + z/3)/2)), where the total mass and mass 
ratio are denoted as rritot = J2p^p = ^p/^tot, respectively. The orbital frequency lo 

for the triangle satisfies cu^ = rritot /o?. 

Henceforth, it is convenient to employ the COM coordinates (X, Y) that can be obtained 
by a translation from (x, y). In the COM coordinates, the location of each mass at any time is 
expressed as Xp = ap{cos{ujt+6p), sm{ujt+6p)), where Op is defined as ai = \/ x\;q^i + Vcom^ 



as = a/ (3i/2a/2 - xcomY + (a/2 - ycouY, and ag = a/x^om + - VcomY, respectively, 
and 9p denotes the angle between the new X-axis and the direction of each mass at t = 
(See Fig. [3]). 
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By using the standard quadrupole formula, the energy loss rate for the Lagrange's orbit 
is expressed as 



dE _ 32 



Vpiyqala^ sin^ (6'p - 6^) 



p<q 



(1) 



The equation of motion for each body is rewritten in an effective one-body form as [2 
dPXp/dt^ = —MpXp/\Xp\^ , where we define the effective mass as 

3/2 



Mp = mi 



ot 



(2) 



The orbital frequency is the same for each body, which provides an identity as Mp/a^ = 
Lo^ from the above effective one-body equation of motion. One can reexpress as = 
{Mp/mtotY^'^a in terms of Mp because tu^ = rritot/o?'- 

For the triangle solution, we obtain the sum of the Newtonian kinetic and potential energy 

as 



m 



tot 



tot 



2a 



2/3' 



(3) 



By assuming adiabatic changes, we use the energy balance between the system energy loss 
and gravitational radiation. We find 



1 da 
a dt 



64 



rritot 



2/3 



'^Ep^q'^p'^q{^, 



2/3 



Mq 

mtot 



2/3 



sm"(6'p - 9g) 



5 a^ 



Y.p^q ^P^q Y.p ^P ( mtlt 



2/3 



-,(4) 



which provides the shrinking rate of the triangle due to gravitational radiation reaction. 

Since the gravitational waves frequency few is twice of the orbital one, we have /q^^ 
rritot/T^'^o?. Therefore, d\n fcw/dt = —{3/2)d\rLa/dt. Using this in Eq. (jl]), we obtain 



1 dfcw 96 



few dt 

where we defined a chirp mass as 

2 



8/3/K.5/3 .8/3 



(5) 



Mchirp — 'nT'tot 



2/3 



mtot 



2/3 



sm" (6'p - 6'g) 



Y^p^q^P^q Y.p^p{^mtot 



2/3 



3/5 



(6) 
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It is worthwhile to mention that the frequency sweep for the triple system can take the same 
form as that for binaries. One can show that Eq. ([6]) recovers the binary chirp mass in the 
limit of m3 0. 

Equation suggests that we cannot distinguish two cases of the binary and triple 
systems by using only the quadrupolar parts even if the frequency sweep is observed. 
Octupole waveforms. — In a wave zone, the gravitational waves denoted by hj^ can be 



37|. The ratio of the octupole part to 



expressed asymptotically in multipolar expansions 
the quadrupole one is of the order of v/c, where f is a typical velocity of the matter. For 
instance, it is about ten percents if a = lOOrritot, which is assumed in order to exaggerate 
the octupole correction in Fig. HI 

After straightforward calculations, one can obtain an expression of octupolar parts of the 
gravitational waves that are generated by the three-body system for the Lagrange's solution 
with arbitrary mass ratio. For instance, one of the relevant octupole moments is expressed 
as 

1 ^ 1 ^ 

Ixxy = -^^f^plXpl^ sm{ujt + 9p) - - ^mp\Xp\^ cos 3{ujt + 9p). (7) 

p=i p=i 
Ixyy can be obtained by interchanges as x ^ y and sin cos. By using such analytic 
expressions, one can obtain the octupole contributions to waveforms. 

It should be noted that no octupole radiation is emitted along the orbital axis for any 
planar motions. Let us take the observational direction along x-axis. Then, we have only + 
mode without x mode. Figure H] shows that a difference between the waveforms (one by the 
binary and the other by the triplet) comes up at the octupole order. The octupole radiation 
amplitude by binaries is proportional to the mass difference jsS]. On the other hand, the 
octupole radiation exists for triangles even if they are all equal masses. Cases of various 
mass ratios and observational directions are a topic of future study. 

Conclusion. — In summary, we have examined different numbers of self-gravitating particles 
in gravitational waves astronomy. In order to track the evolution of the similar waveforms 
from the two-body and three-body systems, we have defined a chirp mass to the three- 
body case. We have shown that the waveforms at the quadrupole level cannot distinguish 
the sources even with observing frequency sweep. Our example suggests that theoretical 
waveforms including higher multipole parts will be important for classification of such similar 
imprints. Higher post-Newtonian corrections both to the waveforms and to the motion of 
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FIG. 4: Gravitational waveforms in arbitrary units for a binary (solid black curve) with mi : 
1712 = 2:3 and a Lagrange solution (dotted red one) with rrii : m2 : ms = 1:2:3, where both the 
quadrupole and octupole parts are included. As a reference, we give the quadrupolar waveforms 
from the same sources (dashed blue). We assume a = lOOmto^ in order to exaggerate a correction 
by the octupole (nearly ten percents expected in this figure). The direction to the observer is along 
X-axis. One can see that the dashed blue curve will overlap with the solid black one after they 
are shifted by choosing the initial phase. This coincidence is because the octupolar waves for the 
binary case are proportional to the mass difference 



and thus relatively small in this figure. 



bodies should be incorporated. This is a topic of future study. In particular, the stability 
of the Lagrange orbit due to general relativistic effects is poorly understood. 

It is conjectured by induction from our result that classification of N (or less) particles 
producing (nearly) the same waveforms requires inclusions of the i-th multipole part with 
i < N. Cases of i < N are realized for instance by the criss-cross and figure-eight. Proving 
(or disproving) this conjecture is left as future work. 
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